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Abstract. We discuss coupling from the past techniques (CFTP) for pertur- 
bations of interacting particle systems on "L d with a finite set of states, within 
the framework of the graphical construction of the dynamics based on Poisson 
processes. We first develop general results for what we call CFTP times with 
ambiguities. These are analogous to classical coupling (from the past) times, 
except that the coupling property holds only provided that some ambiguities 
concerning the stochastic evolution of the system are resolved. If these ambi- 
guities are rare enough on average, CFTP times with ambiguities can be used 
to build actual CFTP times, whose properties can be controlled in terms of 
those of the original CFTP time with ambiguities. We then prove a general 
perturbation result, which can be stated informally as follows. Start with an 
interacting particle system possessing a CFTP time whose definition involves 
the exploration of an exponentially integrable number of points in the graphi- 
cal construction, and which satisfies the positive rates property. Then consider 
a perturbation obtained by adding new transitions to the original dynamics. 
Our result states that, provided that the perturbation is small enough (in the 
sense of small enough rates), the perturbed interacting particle system too 
possesses a CFTP time (with nice properties such as an exponentially decay- 
ing tail). The proof consists in defining a CFTP time with ambiguities for 
the perturbed dynamics, from the CFTP time for the unperturbed dynamics. 
Finally, we discuss examples of particle systems to which this result can be 
applied. Concrete examples include a class of neighbor-dependent nucleotide 
substitution model, and variations of the classical voter model, illustrating the 
ability of our approach to go beyond the case of weakly interacting particle 
systems. 



1. Introduction 

The present paper discusses coupling from the past (CFTP) techniques for inter- 
acting particle systems. The key idea of CFTP, as described in the seminal paper 
[14] by Propp and Wilson, consists in simulating coupled trajectories of a finite 
state-space Markov chain from further and further in the past, until eventually the 
present state of the Markov chain is the same for all trajectories, regardless of their 
starting point. One thus obtains an exact realization of the stationary distribution 
of the corresponding Markov chain, and, under a certain monotonicity condition on 
the transitions of the chain, CFTP leads to a practical algorithm for sampling from 
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the stationary distribution. Many extensions of this scheme have been developed 
since, notably to include processes on more general state-spaces, and situations 
where the monotonicity condition is not met (see the online bibliography [17] ). 

Here, we consider interacting particle systems in the sense of [11] , that is, continuous- 
time Markov processes describing the evolution of a system of states attached to the 
sites of Z d , the evolution at a site being governed by local transition rates involving 
the states of the neighboring sites. Our discussion is limited to particle systems 
with a finite state space whose dynamics can be prescribed by a finite family or 
transition rules (see below for a precise definition) . For an ergodic particle system, 
one is interested in using CFTP to sample from the stationary distribution of the 
system. In general, it is not feasible to sample from the full stationary distribution, 
if only because a full configuration of the particle system is an infinite-dimensional 
object, comprising one definite state for each site of 7L d . A more reasonable goal 
is the following: given a finite set of sites in Z d , use CFTP to sample from the 
marginal of the stationary distribution on this set of sites. This turns out to be 
possible when the interacting particle system possesses what we call in this paper 
a CFTP time, a precise definition being given below. 

In the rest of this introduction, we give a quick overview of our results, some formal 
definitions and statements being postponed to later sections. A discussion of the 
existing literature and how the present work fits into it, is given at the end of the 
introduction. 

1.1. Interacting particle system dynamics and the graphical construction. 

We consider interacting particle systems with a finite state space <S, whose set of 
sites is T, d for some d > 1. To specify the dynamics of the system, we use the notion 
of a transition rule. Such a rule is a triple 9\ = (/, A, r), where A is a finite subset 
of Z d , / : S A S is a map, and r > is a non-negative real number. Given a 
configuration of the system rj = (r](z)) ze ^d G S z , and x G Z d , we denote by fR x rj 
the configuration defined by (%K x rj){x) — f {{rj{x + y)) y& A) and (9\ x r])(z) = rj(z) for 
z/i. (Our convention when A = is that the set S A is a singleton on which / 
takes a single well-defined value.) 

Now, given a finite list of such transition rules (9li)iex = (/i, Aj, r,)j £ i, we consider 
the interacting particle system dynamics specified by the infinitesimal generator C 
satisfying 

(1) C<t>( V ) ^^r, (<fi(m X V) ~ 4>{v)) , 

for all functions <f> : S z — > R that depend only on a finite number of coordinates. 
Informally, |T]) means that, during an infinitesimal time-interval of length dt, in- 
dependently at each site x G Z d , the transformation 5Hf is applied to the current 
system configuration with probability ridt. It is standard to check (see [11]) that 
(fTJ) uniquely characterizes a cadlag continuous-time Markov process (rji)t on <S Z 
equipped with the product topology. 

In the sequel, we assume that the dynamics is in fact built through the graphical 
construction associated with the list of rules {pK%)iex ( see EE] for examples of this 
construction). Specifically, we consider a Poisson point process V on 7L d x I x K 
with intensity J defined by dJ{x, i, t) := ridicjd ®cj® (wjix, i, t), where c%d and 
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ci denote the counting measure respectively on Z d and X, while £r denotes the 
Lebesgue measure on K.. The realization of the point process V prescribes the 
dynamics of the particle system through the fact that, for every x, (r) t (x))t is a 
jump process whose state may change only at times t for which there exists an (a.s. 
unique) i such that (x, i, t) G V, and that, for such a t, one has 

(2) 

Given t\ < ti and £ G S z , define (£) to be the configuration of the system at 
time t<i obtained by starting in configuration £ at time t\— , and using the transitions 
specified by V through ([2]) . We refer to $ as the stochastic flow on S 1 induced by 
V. Given t < 0, we use the notation V t = V n (Z d x [t 3 0[xZ), and let J" t = a(V t ). 

1.2. CFTP times with and without ambiguities. We now consider the cou- 
pling properties of the flow $. We say that a negative and a.s. finite random 
variable T is a CFTP time (for site zero) if the following property holds on the 
even tQ{T>-oo}, 

(3) for all 6,6 G [*°T&)](0) = [*°T (&)](0). 

One checks that the existence of a CFTP time implies ergodicity of the particle 
system. Moreover, starting from an arbitrary configuration £ G 5 Z , the distribu- 
tion of [3>t - (£)](0) * s exactly the marginal at site of the corresponding stationary 
distribution. To obtain a sample from the marginal of this distribution on an ar- 
bitrary finite set of sites, one then has to repeat (after suitable translation) the 
procedure leading to [3>jr~(£)](0) to each site in the set of interest. 

The notion of CFTP time with ambiguities is a weakening of the notion of CFTP 
time, in which property ([3]) holds only when the "ambiguities" associated with 
the rules attached to a specific random subset H of Vt, are resolved. To give a 
precise definition, let us consider, for each a = (x,i,t) G V, and s < t, the random 
variable e(a,£, s) denoting the value at site x produced by the application of the 
rule attached to a when starting in state £ at time s— , more formally: 

(4) e(a, a ,0 = [*&)](*)■ 

When there exist two distinct £1,6 such that e(a, s, £1) 7^ e(a, s, £2)1 w e say that 
there is an ambiguity as to the result of the application of the rule attached to 
a, when we start at time s— . A coupling time with ambiguities consists of a 
negative a.s. finite random variable T, together with a random subset H of Vt, 
which is assumed to be finite on the event {T > —00}, and such that the following 
modification of ([3]) holds: 

(5) for all 6,£ 2 G S z \ [$SJr(6)](0) = [$§T(6)](0) provided that 

e(a,T,£_ 1 ) = e(a,T,6) for all a G H. 

Note that, when H = 0, ([5]) is exactly ©. In addition, we require that H has the 
stopping property, i.e. H nVt is J-j-measurable for all t. 



Throughout the paper, we adopt the convention that all possible pathologies of CFTP times 
and their variants are concentrated on the event that the corresponding time takes the value 
—00, so that one does not have to bother excluding undesirable zero probability events when the 
corresponding time takes a finite value. 
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1.3. Description of the main results. Our first main result is that, starting 
from a CFTP time with ambiguities (T,H), one can build an actual CFTP time 
T* , provided that H contains few enough points on average. To give a precise 
statement, introduce the quantity 



\(x,t,i)£H ) 

Theorem A. If q < 1, one can construct a CFTP time T* for the interacting 
particle system. 

The construction of T* is explained in Section [2j Here is an informal description. 
Starting with Ambo := {(0,0)}, we recursively define a sequence (Amb„)„>o of 
random subsets of Z d x R in the following way. First, we apply the coupling time 
with ambiguities (T,H) at each space-time point in Amb„. This generates a set 
of elements of V , with respect to which ambiguities have to be resolved. Then 
Amb„ + i is defined as the set of space-time points upon which the resolution of 
these ambiguities directly depends, i.e. for a — (x,i^t), the set {(x + y,t); y £ A;}. 
The overall set of points generated by this process is Amboo :— lj n>0 Amb„, and 
T* is defined as the lowest value of T obtained when applying the coupling time 
with ambiguities (T, H) to the space-time points in Amboo ■ The idea underlying 
this construction is that, if Amboo is finite, one can resolve ambiguities in a step- 
by-step manner, starting from the points in Amboo that are furthest in the past and 
thus associated with an empty set of ambiguities, down to the origin where we can 
determine the value of (£)](0) ( a precise formulation is given in Proposition l5.ll 
in Section [5]). 

Without giving precise statements (see Theorems ICl and iDl in Section [3]), let us 
mention that, in addition to Theorem lAl it is possible to obtain estimates on the tail 
of T* and on the range of its space-dependence (in terms of bounds on exponential 
moments), from analogous properties for T and H. 

Our second main result deals with perturbations of interacting particle systems. To 
formalize this notion, consider a particle system whose dynamics is defined by a list 
of rules (9ii)iei"' This corresponds to the original, unperturbed, particle system. 
Then consider the dynamics defined by a list of rules of the form ($Hi)igX"uiPj 
where I p is disjoint from X u . This corresponds to the perturbed particle system. 
Our result gives general conditions under which the existence of a CFTP time T u 
for the unperturbed particle system leads to the existence of a CFTP time for the 
perturbed particle system, provided that the perturbation is small enough. Our first 
condition is that the unperturbed dynamics possesses the positive rates property, 
which means that, for every v £ S, there exists a rule with index in I u whose 
application inconditionally leads to the value v. Our second condition requires that 
the definition of T u involves the exploration of an exponentially integrable0 number 
of points in V, a notion whose precise formulation is given in Section[2j and involves 
what we call the exploration process associated with T u . Finally, the smallness of 
the perturbation is measured through two parameters e and k, that admit explicit 
definitions in terms of (9ti)i62> an d (9^i)iexp (see Section [2T4|) . 

2 For a non- negative random variable X, we say that X is exponentially integrable if there exists 
> such that E{e flX ) < +oo. 
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Theorem B. Assume that the unperturbed dynamics has the positive rates prop- 
erty, and possesses a CFTP time associated with an exploration process whose total 
number of points is exponentially integrable. Then, for any perturbation with small 
enough e and k, we can construct a CFTP time with ambiguities (T, H) for the per- 
turbed dynamics, satisfying condition ([6]), so that the corresponding T* is a CFTP 
time for the perturbed dynamics. 

The construction of (T, H) is explained in Section[3] Note that, in addition to satis- 
fying (|6]), (T,H) also satisfies the assumptions of Theorems [Cl and [Dl for sufficiently 
small e and k, leading to exponential moment bounds on the tail of T* and on the 
range of its space-dependence. This extension of Theorem [B] is stated as Theorem 
Ein Section El 

We illustrate Theorem [B] with applications to several kinds of interacting particle 
systems. A first class of examples is given by systems that satisfy what we call 
the finite factor property (see Section |4}. Informally, this means that the state 
of a site at a certain time depends on the points in the graphical construction 
and on the initial condition only through a window of fixed size around x. The 
simplest example is provided by dynamics for which distinct sites do not interact, 
whose perturbations correspond to weakly interacting interacting particle systems. 
A more elaborate concrete example is a family of nucleotide substitution models 
called RN+YpR (see [T]), which allows for arbitrarily strong interactions between 
neighboring sites and yet satisfies the finite factor property. Another class of sys- 
tems to which we apply Theorem [B] (and for which the finite factor property is not 
satisfied), is given by noisy voter models on 7L d . Specifically, we consider the case 
of a classical linear voter model with an arbitrary finite alphabet, and a variation 
we call the voter model with asymmetric polling, that uses the classical {+, — } 
alphabet (see Section 2]). Note that these examples too extend beyond the weakly 
interacting case. 

Although we do not enter into the details here, let us note that the existence of a 
CFTP algorithm is not only interesting for simulation purposes, but can also pro- 
vide interesting theoretical results on the particle system. Indeed, the existence of 
a CFTP time automatically implies that the interacting particle system is ergodic, 
and estimates on the tail of the CFTP time such as those provided by Theorem 
[Cl immediately lead to bounds on the speed of convergence to the stationary dis- 
tribution. Similarly, estimates on the range of the space-dependence such as those 
provided by Theorem [D] yield bounds on the decay of spatial correlations. As a 
consequence, our results can be readily used to derive interesting conclusions about 
the behavior of the perturbed particle systems to which Theorem [Bl applies. 

Also, note that, in this paper, we do not explicitly address issues related to the 
practical implementation of CFTP. However, from the definition of T* in terms of 
T and H it should be clear that, if T and H lend themselves to an actual algorithmic 
implementation, so is the case for T* . Similarly, Theorem [Bl yields an actual CFTP 
algorithm for the perturbed particle system provided that T u and the associated 
exploration process are compatible with an actual algorithmic implementation. 

Finally, let us point out that a key role in the proof of our results is played by first- 
moment arguments, that allow us to essentially bypass the quite intricate analysis of 
dependencies that would otherwise be required to study the combinatorial objects 
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we have to deal with (e.g. the sequence of sets (Amb n ) n >o, or the tree-indexed 
exploration process X used to define a coupling time with ambiguities from T u ). 
In fact, first-moment arguments allow us to largely ignore these dependencies and 
obtain results in very much the same way as for classical branching processes. 

1.4. Discussion. For ergodic particle systems satisfying a monotonicity condition 
similar to that of [14], CFTP is always possible, as shown by van den Berg and 
Steif in [16] . For systems lacking monotonicity, CFTP algorithms have been devel- 
oped under "high- noise" or "weak interaction" type assumptions, meaning that the 
strength of the interaction between neighboring sites has to be sufficiently small. In 
other words, the particle system under consideration has to be a sufficiently small 
perturbation of a system in which distinct sites do not interact. One example is 
given by Haggstrom and Steif [10] , who use a bounding set approach to control the 
coalescence of trajectories (see also [3] for some refinements). Another example is 
given by Calves, Garcia, Locherbach [7] (see also [E1IH]), whose approach is based 
on a branching constructior0 of which the one we use in the present paper can be 
seen as a generalization. One interest of the present paper is that it provides a 
general criterion under which small perturbations of an interacting particle system 
retain some of the CFTP properties of the original unperturbed system, allowing 
one to go beyond the weakly interacting case. Let us mention that some of our 
applications overlap with the recent paper |13j , where a specific kind of perturba- 
tion of noisy voter models is considered, and ergodicity is proved for sufficiently 
small perturbations. Finally, note that the present paper is a revised and extended 
version of an earlier manuscript [2] , where the results were limited to perturbations 
of RN+YpR nucleotide substitution models. 

1.5. Organization of the paper. The rest of the paper is organized as follows. 
Section [2] contains the definitions of the various notions and objects that were 
encountered in the introduction but not formally defined. Section [3] contains ad- 
ditional results that were not stated in the introduction. Section 2] contains the 
examples of application of Theorem [B] Section [5] describes the proofs of Theorems 
lAllCl andlDl Section[5]describes the proof of TheoremlBl (and its extension Theorem 

H). 

2. Some formal definitions 

2.1. Canonical probability space. We assume throughout the paper that V is 
defined on a canonical probability space (f^J 7 , P) that we now describe. First, 
is the set of locally finite subsets ui of 7L d x 2" x M. satisfying the additional 
requirements that (i) no two points in u) share the same M— coordinate, (ii) for 
every (x, i) S Z d x I, both sets lo n ({(x,i)} x M + ) and uj n ({(x,i)} x R_) are 
infinite, (iii) for every (x, i, t) 6 u, and any sequence (y n ,jn, s n )n>o starting at 
(ya,ja,s ) = (x,i,t) and satisfying y n+1 <E y n + A jn and s n+ i < s n for all n > 0, 
one has lim„_i. +00 s n = —oo. Then J- is defined as the cr— algebra generated by all 
the maps of the form cj H> | u fl B \ , where B is a Borel subset of Z d x I x R. Finally, 



very similar construction was already used in [6] to devise CFTP algorithms in a different 
framework. In fact, various constructions of this kind appear in the literature, though not explicitly 
in the context of CFTP, see e.g. [5j [4] . 
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we set V(lu) := w, and P is uniquely denned on (17, J 7 ) by the requirement that V 
is a Poisson process with intensity J (this definition makes sense since (i) , (ii) , (iii) 
are almost sure properties of a Poisson process with intensity J) . Note that thanks 
to (i), (ii), (iii), 3>t°(£) is well-defined for any t\ < t%, £ G S z and w 6 17. 

Given € Z rf x R, we define the space-time shift r^t on 17 by r X) t(a;) := 

U( z i s ) a {( z - s — t)}. (This definition is possible since (i), (ii), (iii) all are shift- 
invariant properties.) Given a G NU{±oo}, we define V~ a ' a := Vn([—a, a] d xlxM.), 
and let T~ a ' a = a{V^ a ^ a ). On occasions, we use the notation Vf a,a := V n 
{[-a,a] d x [t,0[xl). 

Finally, to properly define the notion of a random subset of V, we introduce the 
space 17 formed by the subsets of elements of 17, equipped with the a— algebra 
J- defined just as T . Note that, as a rule, in the rest of the paper, we mention 
measurability issues only when they involve a non-trivial argument. 

In the sequel, we have to consider two distinct probability spaces, one associated 
with the unperturbed dynamics, and one with the perturbed dynamics. We denote 
by (17 u , J 7 ' 1 , P u ) the probability space associated with the unperturbed dynamics 
specified by (jHj)iexu, while (17, J 7 , P) refers to the perturbed dynamics specified by 
the full list (Pli)iGZ. The corresponding Poisson processes are denoted respectively 
V u and V. 

2.2. Construction of T*. Given a coupling time with ambiguities (T,H), define 
by induction the following random sequence of subsets of Z d x] — oo, 0]: 

(7) Amb :={((), 0)}, Amb„ +1 := (J (J \J{(z + y,s)}. 

(ii)£Amb„ ( z i s ) gT -i offoTx t yeAi 

Then let Amboo := \J n>0 Amb„, and T* :— in^x^eAmb^ t + T o r x> t in the case 
where Amboo < +oo, while T* := — oo otherwise. 

2.3. Exploration process. Here, we define the notion of an exploration process 
attached to an interacting particle system. 

Given a non-empty finite subset B C Z d and t < 0, define n(B, t) to be the element 
of PC\(B xlx] — oo, t[) with the highest time coordinate (this is always well-defined 
with our choice of 51). Let fi/ denote the set of all finite subsets of sets in 17, and 
let 9 denote a measurable map from 17/ to the set of finite subsets of Z d (where 17 f 
is equipped with a a— algebra defined as T). 

The exploration process associated to 9 is defined as follows. Start with Xo := 0, 
70 := 0. Then, for all n > 0, let B n := 9(X n ). If B n ^ 0, denote n(B n ,j n ) =: 
(x n ,i n ,t n ), and let X n+ i := £„ U {(x n , i n , t n )} and 7„+i := t n . If B n = 0, then 
X„+i := X„,7„ + i := j n . Let Xoo := \J n > X n and 700 := lim n _>. +00 j n . The total 
number of points in the exploration process is then defined as |Xoo|. 

We say that a CFTP time T defined on (17, J 7 , P) is associated with such an ex- 
ploration process if, on {|Xoo| < +00}, one has T — 700, or equivalently, T = 
inf{i; (x,i,t) £ I m }, while T = —00 when |Xoo| = +00, and if, on {T > —00}, 



Note that we are dealing with negative numbers here, so that the highest time coordinate here 
corresponds to the time coordinate with the least absolute value. 
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the value of (£)](0)> which is the same for every £ £ , is measurable with 
respect to Xoo. 

We shall always assume that there exists a deterministic function f3 : N — > N such 
that for every I > 0, 

(8) 0(x t ) c {-m,---,m} d , 

and such that (5(1) — 0(£) as £ goes to infinity. 

In the sequel, we assume that a CFTP time T" for the unperturbed dynamics is 
defined on (fl u , T u , V u ), and that T u is associated with an exploration process of 
the type we have just described. We denote the corresponding process by X u to 
emphasize the fact that this process is defined for the unperturbed dynamics, but, 
for the sake of readability, we use 9, x n , i n , t n , j n , etc. instead of the awkward 9 U , 

Remark 2.1. Given a CFTP time T and an exploration process (X n )n> the fact 
that T — 7^ does not in general imply that [$XT(£)](0) is measurable with respect 
to Xoo, so this last condition has to be added to the definition of an exploration 
process associated to a CFTP time. 

2.4. Positive rates property, e and k. The positive rates property for the set 
of non-perturbative rules (^Hi)igi" means that, for every v € 5, there exists a rule 
with index in I u which is of the form (A, /, r) with r > 0, A = and f = v. 
We denote by l v theQ index of such a rule. We then control the smallness of the 
perturbation of ($Hi)i£i u by (^Hi)iei p through the following two parameters 



(9) e := sup 

v£S 



(10) 



,je2>; vefj(A 




2.5. Construction of a coupling time with ambiguities (T, H) from T u . We 

now define what we call the exploration process with locking of perturbative ambi- 
guities attached to the perturbed dynamics, and associated to the map used to 
define the exploration process X u of the unperturbed dynamics. This is the process 
we use to define a coupling time with ambiguities for the perturbed dynamics. 

Informally, the construction can be described as follows: run the exploration process 
associated with the unperturbed dynamics on V (instead of V u ). When an a = 
(x,i,t) corresponding to a perturbative rule, i.e. i € I p is encountered, split the 
exploration process into exploration processes evolving in parallel, one for 

each v 6 f(Ai), in which (x, i, t) is replaced by (x, i v , t). 

The formal construction uses a recursively constructed tree T to label the process. 
Let us start with the root of T, denoted r, for which we set X r := and 7 r := 
0. Then, for any vertex a £ T, we set B a := 9(X a ). Then, if B a ^ 0, denote 
n(-B a ,7 a ) = (x a , i a ,t a ). If i a £ I", we attach a single son b to a, and let Xb := 



^One may assume without loss of generality that, for any v, there is a single such rule, since 
identical rules with distinct indices may always be merged into a single rule. 
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X a U {(x a ,ia,ta)} and 76 := t a . On the other hand, if i a G I p , we attach to a a 
list of sons (&„, v G fi a {S Aia )) , and let Xf, v := X a U {(^a, &D, i a )} and 7^ := t a . If 
£? a = 0, then a has no son. Finally, we let Xoo := UaeT 
We now define T by 

(11) T := inf{£; (x,i,t) G Xoo} if |Xoo| < +00, 
while T := — 00 if |3£oo| = and H is defined by 

(12) H := {(x a ,i a ,t a ); a G T', 

where T' denotes the subtree of T obtained by removing the leaves of T. Note that 
one may view X^ as the result of the exploration process associated to a certain 
map 9. However, the representation with a labelled tree turns out to be more 
convenient for our purposes. 

Remark 2.2. The definition of X^ makes sense whether or not \3t% \ has an 
exponentially decaying tail, as assumed in the statement of Theorem\B^ provided that 
\Xoo\ is a.s. finite, (T,H) is indeed a CFTP time with ambiguities (see Proposition 
\6.1\) . and the role of the tail condition on |3£^,| in Theorem\B\is indeed to ensure 
that \Xoo\ is a.s. finite. 

3. Additional results 

The first result shows that the exponential moments of T* can be controlled in 
terms of similar moments for T and H . 

For A G R, define 

(13) Ay (A) := E(exp(AT)), A H ,u m e(\) := E f £ \A\ cxp(At) J . 

Theorem C. Assume that A is such that A^(A) < +00 and Ajj )t j me (A) < 1. 
Then 

E(exp(AT*)) < Ar(A)(l - A J? , time (A))- 1 . 

Our next result deals with the range of space-dependence of T* . To formalize this 
notion, say that a NU {+oo}-valued random variable L defines a stopping box 
in Z d if, for any a G N, one has {L = a} G T~ a,a . We say that an a.s. finite 
such random variable bounds the width of a CFTP time T if, on {T > — 00}, 
the value of [3>jT (£)](Q) (which by definition does not depend on the choice of 
£ G S zd ) is measurable with respect to T~ L ' L . We say that L bounds the width 
of a CFTP time with ambiguities if H is measurable with respect to T~ L ' L and 
if there exists a measurable map3 O such that, on {T > —00}, for all £ G <S Z , 
[*t" (0] (0) = 6 (L, V- L > L , (e(a, T, £)) aeJ? ) . 



"To be more specific about measurability assumptions concerning (e(a, T, , we assume that 

it is encoded as the random subset of Z d X I X R X S defined by U a gjj (<*> e ( Q i ^1 ?))> where the 
a— algebra on the set of locally finite subsets of lA X X X R X S is generated by maps of the form 
ro — > \w n S|, where i? is a Borel subset of Z d X X X R X S. 
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For A € R, and 1 < q < d, define 

(14) A L (A) := E(e AL ), A H „(X, q) := E I ]T ]T e M*,+«.) j . 

Then define L^ by L^ := sup^^ sup (2 , t)eAmboo x q + L o r x ,t, and LI by LI := 
infi< g < d inf (;r . it ) eAmboo x q — Lo r Xjt . Finally, let L* := max(L+, -L~). 

Theorem D. If Q < 1, and L bounds the width of (L, H) then L* bounds the width 
of T* . Moreover, if A > is such that A^(A) < +oo and A# )Space (±A, q) < 1 for 
all q, then 

E(exp(AL;)) < A L (A) sup (1 - A H ,space 

l<q<d 

E(exp(-AL*_)) < A L (A) sup (1 - A H , space (-A, q))~ l . 

l<q<d 

Finally, we have the following extension of Theorem [5] 

Theorem E (Extension of Theorem [Bjl . Under the assumptions of Theorem\B\ for 

any list of perturbative rules with small enough e and k, the pair (T, H) defined by 
(|lip and (|12p defines a CFTP time with ambiguities that satisfies the assumptions 
of Theorem\^ together with the assumptions of Theorems\Qand\J^for small enough 
\X\. 

For the sake of readability, we did not include explicit estimates in the statement of 
Theorem [El However, looking at the proofs given in Section [6j it is easy to obtain 
explicit control upon the characteristics of {T,H) (namely, g, At, A#, time , A^, 
Ajy l5p ace) in terms of e, n, and the parameters (^Ri)ii£X u of the unperturbed model. 
Combined with Theorems [5] O and [12 this leads to an explicit control upon the 
characteristics of the resulting CFTP time T*. 

4. Applications 

In this section, we give some examples of dynamics which satisfy the properties re- 
quired for the unperturbed dynamics in Theorem [Bl Since we discuss unperturbed 
dynamics only, it is unnecessary to use " superscripts to distinguish between per- 
turbed and unperturbed dynamics, and consequently such superscripts are not used 
in this section. 

4.1. Perturbations of finite factor models. We say that the dynamics possess 
the finite factor property if there exists b £ N such that, for all t < 0, [3 > t _ (0] is 
measurable with respect to ~p~ b ' b and (£,(x); x e {— b, . . . , 

Proposition 4.1. Any dynamics with the finite factor property and the positive 
rates property satisfies the assumptions of Theorem\B^ i.e. there exists a CFTP time 
associated with an exploration process whose total size has some finite exponential 
moment. 

Proof. Given X € Qf, let q :— \X\, and denote by (yk, jk, Sk)o<k<q-i the list of 
elements of X, indexed by decreasing order of time, so that sq > ■ ■ ■ > s 9 _i. Let 
also h := (2b + l) d . Now consider the exploration process associated with the map 
defined as follows. Set 9(X) :— when the following three conditions are met 
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a) \X\ > h, 

b) {y q -i,...,y q - h } = {-&,..., b} d , 

c) for all q - h < k < q - 1, A jk = 0, 

Otherwise, let 6(X) := {— b, ...,b} d . Denote by (X n ) n the corresponding explo- 
ration process, and observe that condition ([8]) is satisfied with f3(£) :— b for all 
I. One checks that given X n , the probability that % n +h satisfies 8(X n +h) = 
is bounded below by h\h~ h p h , where p := (Eiez 1 "* 1 ^ = 0))(Eiez r <) _1 - This 
proves the fact that there exists /i > such that E(e M ' Xoo |) < +oo. Let us now 
check that T := 7^ is a CFTP time for the dynamics, associated with the explo- 
ration process defined by 9. Define U to be the a.s. finite smallest index k such 
that 8(Xk) — 0- From conditions a) b) c), one has that, on {U < +00}, for all 
x G {— 6, . . . , b} d , [$y t/_h (£)](x) takes the same value for every £, and this value 
is measurable with respect to (yu,ju), (yu-h+i, ju-h+i)- On the other hand, 



-6,6 



lu- 



the fact that 9(Xk) — {—b, . . . , b} d for all k < U — h shows that Xu-h = V. 
Now by the definition of the flow, one has that 

[*t"(0](o) = [*£_ fc (*F- h -(Q)] (o). 

The finite factor property then yields that [^ _ (^)](0) is the same whatever the 
value of £, and that this value is measurable with respect to Xjj- Q 



The simplest example of dynamics with the finite factor property is the case where 
distinct sites do not interact, i.e. Ai C {0} for every i £ 2. In this case, the 
factor property holds with 6 = 0, and each site evolves independently according to 
a continuous-time Markov chain on S. 

A more sophisticated example, whose study was our original motivation for this 
work, is the so-called class of RN+YpR nucleotide substitution models, see PQ, 
whose goal is to provide tractable models that include neighbor-dependent effects 
such as the well-known hypermutability of CpG dinucleotides. These models use 
the nucleotidic alphabet S := {A, C, G, T} as their state space, and Z as their set of 
sites, with S z being an idealized representation of a DNA sequence. Additionally, 
S is divided into the set of pyrimidines Y := {C, T}, and purines R := {A, G}, and 
we say that Y is the type of C and T, while R is the type of A and G. 

The RN+YpR dynamics is specified through the following list of rules (each rule is 
of the form (/, A, r)): 

• unconditional rules: for each v G 5, a rule with A := 0, r > and f = v; 

• transversion rules: for each v £ S, a rule with A :— {0} and f(w) := v if v 
and w are not of the same type, f(w) :— w otherwise; 

• transition rules: for each v G S, a rule with A := {0} and f(w) := v if v 
and w are of the same type, f(w) := w otherwise; 

• left-dependent rules: for each u G Y, v G R, v' G R, a rule with A := 
{-1,0}, f(w-i,w ) ■= v' if (w-i, w ) = (u,v), f(w-x,w ) := w otherwise; 

• right-dependent rules: for each u G Y, v G R, u' G Y, a rule with A := 
{0,1}, f(w , Wi) :=u' if (w ,Wi) = (u,v), f(w ,wi) := w otherwise. 

It turns out (see [I]) that RN+YpR models have the finite factor property with 
b := 1. Let us insist that the rates of left- and right-dependent rules, whence 
the strength of the interaction between sites, may be arbitrarily large, so that the 
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RN+YpR class contains models that are not weakly dependent. Note that one 
can generalize this class of models to produce interacting particle systems with an 
arbitrarily long range of dependence, where the minimal b for which the finite factor 
property holds can be made arbitrarily large, although these seem less biologically 
motivated. Note also that, in the case of the RN+YpR model, one can define 
alternative coupling times which, as opposed to the one defined in the proof of 
Proposition 14. 1\ do not get larger and larger when the interaction strength (given 
by the rates of the rules involving interactions between neighboring sites) gets large, 
see 0. 

4.2. Perturbations of voter-like models. Wc now describe how Theorem [Bj 
can be applied to variants of classical interacting particle systems such as the voter 
model on Z d (see [12]). 

4.2.1. Classical linear voter model. Let p(-) denote a probability measure on Z d 
with finite support. The dynamics of the classical voter model can be defined 
thanks to the following set of rules: 

• state-copying rules: for each x in the support of p(-), a rule with A := {0, x} 
and r :=p(x), with f(wo,w x ) :— w x . 

One might interpret this model as describing the evolution of opinions of individuals 
attached to the sites of Z d , with S representing the set of possible opinions. The 
individual at x waits for a unit exponential time, then chooses a random location 
y £ 1 d with probability p(y — x), and adopts the opinion of the individual attached 
to site y. As such, the voter model does not satisfy the assumptions of Theorem 
IbI since it does not enjoy the positive rates property. As a consequence, we add to 
this model a list of unconditional rules so as to enforce this property: 

• unconditional rules: for each v £ S, a rule with A := 0, r > and / = v. 

We call the resulting model " noisy voter model" . Note that this addition dramati- 
cally changes the dynamics of the voter model, since it automatically turns it into 
an ergodic interacting particle system. Note that we may consider this addition as 
part of the perturbation of the original voter model we want to study, but this part 
of the perturbation has to be included in the dynamics prior to the application of 
Theorem [Bj 

Proposition 4.2. The noisy voter model satisfies the assumptions of Theorem W[ 

Proof. The corresponding exploration process is defined as follows. First 0(0) := 
{0}. Then, given a non-empty X € fi/, denote by (y, j, s) the element of X with 
the lowest time-coordinate. Then let 6(X) := if Aj = 0. Otherwise, Aj is of the 
form {0,x}, and we let 6{X) := {y + x}. We denote by (3£ n ) n the corresponding 
exploration process. Note that condition © is satisfied with (3(1) := sup{|z|; p(z) ^ 
0} x I. To prove that the |£oo| has some finite exponential moment, note that, 
conditional upon X„, if 0(X n ) ^ 0, the probability that the next point to be included 
in X n +i corresponds to a rule of the form / = v for some v, is bounded below by the 
ratio (X^ues r ^ )(Tl l iei r i) 1 - Since in this case 6*(X„+i) = 0, a geometric upper 
bound holds for the tail of |£oo|- Let us now check that T defined as the least 
time-coordinate of an element in Xoo is a CFTP time for the dynamics, associated 
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with the exploration process defined by 9. Indeed, it is clear from the definition of 
the dynamics that if the element of Xoo with the least time-coordinate is associated 
with the rule i v , then (0) = v for all £. □ 

4.2.2. Voter model with asymmetric polling. We now consider a variation upon the 
classical voter model. Let . . . ,A^ m ^ denote a list of finite non-empty subsets 
of Z d , rW, . . . , r^™) denote a list of non-negative real numbers, and take as a state 
space S := {+,—}. The set of rules characterizing our model is the following: 

• polling rules: for each 1 < i < to, a rule with A :— A™ and f(w) := + if 
w x = + for at least one x £ A^\ f(w) := — otherwise. 

We call this model the voter model with asymmetric polling. Here, an individual 
performs a poll within a randomly chosen finite subset of individuals, and adopts 
an opinion that depends on the results of the poll in an asymmetric way: indeed, 
the individual will adopt the opinion denoted 4- if any of the individuals in the poll 
expresses the opinion +, while, to adopt the opinion denoted — , consensus within 
the poll is required. 

As in the case of the classical linear voter model, we add to the above set of rules 
a list of unconditional rules ensuring the positive rates property: 

• unconditional rules: for v — +, — , a rule with A :— 0, r > and / = v. 

The resulting model is called noisy voter model with asymmetric polling. 

Proposition 4.3. The noisy voter model with asymmetric polling satisfies the as- 
sumptions of Theorem ["51 

Proof. The corresponding exploration process is defined as follows. First 0(0) := 
{0}. Then, given a non-empty set X £ 0/, denote by (y,j,s) the element of X 
with the least time-coordinate. Then let 6{X) := if Aj = and fj = +. If Aj = 
and fj = — , then let 9(X) := 9(X \ {(y,j, s)}) \ {j}. Otherwise, Aj is of the form 
A( fc ) for some 1 < k < to, and we let 9(X) := 9{X \ {(y,j,s)}) U (y + A^). We 
denote by (X n ) n the corresponding exploration process. Note that condition ([8]) is 
satisfied with f3(£) := sup{|z|; z e ^\<k<mA^} x t 

To prove that the number of points in \Xoo\ has some finite exponential moment, 
note that, conditional upon £„, if 9(X n ) ^ 0, the probability that the next point 
to be included in corresponds to the rule with / = +, is bounded below 

by the ratio r L+ (X^gz r i) ■ Since in this case 9{X n+ i) — 0, a geometric upper 
bound holds for the tail of |Xoo|. Let us now check that T defined as the least 
time-coordinate of an element in is a CFTP time for the dynamics, associated 
with the exploration process defined by 9. Indeed, it is clear from the definition of 
the dynamics that if the element of Xoo with the least time coordinate is associated 
with the rule t+, [$jT (£)] (0) = + for all £, while, if this element is associated with 
the rule t_, [$§T(f)] (0) = - for all f. □ 

5. Proofs of Theorems [Al ICl and [Pi 

We start with a proposition showing that, if T* is finite with probability one, then 
T* is indeed a CFTP time. 
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Proposition 5.1. 7/P(T* > -oo) = 1, then T* is a CFTP time. 



Proof. Note that, with our definitions, T* > — oo implies that |Amboo| < +oo. 
The proof is by induction on |Amboo|. Specifically, we shall show for all n > 
that the following property (P n ) is true: on HAmbool = n, T* > — oo}, for all 
£i,6 G S % \ [$t*(Ci)](0) = [$^(6)](0). Assume throughout that T* > -oo, and 
let us start with n = 1. If |Amboo| = 1, a first possibility is that H = 0. In this 
case, the definition of a CFTP time with ambiguities shows that [<E>^ — (^i )] (0) = 
[*t"(&)](0) for all & ,fa € S z , whence, since by definition T* < T, the fact that 
[*t"(Ci)](0) = [$°-(£ 2 )](0) for all Ci, 6 €SP*.1fH^ 0, the fact that lAmb^l = 1 
shows that every a — (x, i,t) € H is such that A; = 0. In this case, for any s <t, 
e(a,£,s) depends neither on £ nor s, so that again [$y7(£i)](0) = [ < f>!^(£2)](0) for 
all £i,^2 € 5 Z . We now show that (P n +i) is valid provided that (Pk) is valid for 
all 1 < k < n. Assume that |Amboo| = n + 1. It is enough to prove that, for any 
a = (z,i,t) G such that Aj 7^ 0, e(a, T*,£) admits the same value for every 
£ G iS z . Consider such an a = (z, i, t), let y G Aj, and x := z + y. Then observe 
that, by definition, 

(15) t~1 o Amboo o t x i C [J Ambfc. 

k>l 

Since Amb = {(0,0)} while (0,0) ^ lj fe>1 Ambfc, our assumption that lAmb^l = 
n+1 implies that | Ufc>i Ambfc | = n. We thus deduce from (fT5|) that lAmbooor^ tl < 
n. Moreover, (fT5|) shows that t + T* o r x ,t > I 1 *, so that our assumption that 
T* > —00 implies that T* o r x .t > —00. Our induction hypothesis then implies 
that for all £i,£ 2 G S 2 *, (0) o Tx , t = [$t7(6)](0) o r x . t , which rewrites as 

[*|; T . OT . it Ki)](a:) = [*£t«ot.,« (&)](*)■ We have seen that T* < t + T* o Tx , t , so 
we can deduce that e(a,T*,£) does not depend on £ G <S Z . □ 

Let M denote the intensity measure of the set Ambi, i.e. the positive measure on 
Z d x R defined for all Borel set C by 

M(C) := E(|Ambi n C|) = E ^ ^ l((z + y, s) G C) 

\(z,i,s)£H y€Ai 

We use the notation * for the convolution product of measures on Z d x R. For all 
n > 0, M* n denotes the convolution product M * • ■ • * M with n factors, with the 
convention M*° := <S(o,o)- 

Our key first-moment estimates are given in the next proposition and its corollary. 
Proposition 5.2. For any measurable f : Z d xl4 R +; and any n>0, one has 



E no) < / 



f(()dM™(0, 



with the convention M*° := <Wo)» 
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Corollary 5.3. For measurable f : Z d xR — > R_|_, any non-negative J-q — measurable 
random variable D, and any n > 0, one has 



e( Yl KQ'DotA <e(d)- ff^uixr"^). 



with the convention M*° := <5(o,o)- 

The proof makes use of the so-called refined Campbell theorem (see [15]), which we 
quote here in the special form we need: 

Theorem F. For any measurable map : (Z d xlxl) x f2 — » R + , one has the 
following identity: 



eI ^*(a,P) j = J E(tt( 



>(a,VU {a}))dJ(a), 
where J denotes the intensity measure ofV. 

A crucial property of the sets Amb n is that they possess the stopping property, as 
stated in the following lemma. 

Lemma 5.4. For all n > 0, Amb„ has the stopping property, i.e. for all t < 0, 
Amb„ n Vt is J- immeasurable. 

Before we prove Lemma 15.41 we need the following lemma. 

Lemma 5.5. There exists a measurable map S from Qx] — oo,0[ to f2 such that, 
for all t < 0, 

Hnv t = z(t,v t ). 

Proof. We first prove that there exists a measurable map V from f2x] — oo, 0[ to Q 
such that, for all t > 0, 

(16) HnV t = V(t,V). 

For each (x,i) e Z d x I, let (ip(x, i, fc))fc>i denote the successive points V whose 
coordinate on 7L d x X is (x,i), in decreasing order of the R— coordinate. We let 
ip(x, i, k) =: (x, i, s(x, i, k)). Given a Borel set B of 1 d x I x R, one has that, for 
all t < 0, 

\H n V t n B\ = Y m>(x, i, k) e H)l(s(x, i, k) > t)l(i(}(x, i, k) e B). 

Since H is a measurable map from (Vt, J 7 ) to (VI, F), one has that uj h-> l(ip(x, i, k) S 
H)(ui) is a measurable map from (VI, J 7 ) to R. This is also the case for uj H» 
l(ip(x,i, k) € B)(lu). Finally, (t,oj) >-> (s(x,i,k)—t) is measurable from fix] —00, 0[ 
to R, so this is also the case for l(s(x, i, k) > t). We conclude that (ui, t) i—^HnVt 
is measurable from fix] — 00, 0[ to VI, whence the existence of V. Consider now an 
arbitrarily fixed element ljq G Vt that contains no point with R— coordinate, and 
define the map a from fix] — oo, 0[ to f2 by 

a(u, t) := (0 n (Z d x I x [*, 0[)) U r ,t(w ). 
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One checks that a is measurable by writing 

\a(u,t)r\B\ = ^l(s(a:,*,fc)(w)>t)l(^(x,i,A:)(w)GB) + |-nj,i(£Jo)nB| J 

x,i,k 

where we have extended the definition of ip(x,i,k) to f2 in the obvious way, with 
the convention that s(x,i,k) takes the value — oo when the value of k excesses 
the number of points to be indexed. Now, since, for any given t < 0, H n Vt is 
Tt~ measurable by assumption, one has that, in view of (|16l) . for any t < 0, 

V(t,V) = V(t,a(V t ,t)). 

As a consequence, the conclusion of the proposition is achieved by defining 

E(t,u) :=V{t,a(u>,t)). 

□ 

Proof of Lemma \5.4\ We re-use the notations introduced in the proof of Lemma 
15.51 The proof is by induction. For n = 0, the result is obvious since Ambo := 
{(0,0)}. Now assume the result to be true for a given n > 0. Define D n := 
{(w, i, k); n(if)(w, i, k)) n Amb„ ^ 0}. By definition, one has 

(17) Amb n+1 := (J it (^(x.i,*) ° H ° r x, S (x,i,fe)) > 

(i,i,f:)eD„ 

with the slight abuse of notation that, given a subset C of Z d x I x R, 7r(C) := 

Now consider i < 0. From Lemma 15. 5i we deduce that the map on £1 x [t,0] 
defined by (ui, s) i-> (t^J o Jf o t XjS (u})) R 'Pt(w) is ® 0])— measurable. On 
the other hand, our induction hypothesis shows that for any (x, i, k), the event 
{(x, i, k) £ s(x,i,k) > t} is J^— measurable. We can then deduce from (fl7|) 
that Amb„ + i n "P t is J^— measurable. □ 

Proof of Provosition \5.2[ We proceed by induction. For n = 0, the result is imme- 
diate since by definition Ambo : = {(0,0)}, while, for n = 1, the result is a direct 
consequence of M being the intensity measure of Ambi. Now consider n > 1, and 
note that, by definition, 

(is) E /(O < E 0*,*> 

CGAmb„ + i (i,t)eAmb„ 

with 

9x,t ■= E E /( z + 2/'*' s )- 

(a,»,s)€T~JoJ?OT !1 . it y£-4i 

For (a;, z, t) £ Z d x 2" x M, define 7r(x, i, t) := \J zE a i( x + z ' *)}■ Now we rewrite 

E = E 

(x,t)eAmb„ a£P 

with 

*(a,w) := l(?r(a) C Amb n (u;)) E 0*,tM- 

(x,t)£7r(a!) 
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Applying Campbell's theorem (Theorem IF]) . we deduce that 
(19) E 



E /(C)) < Jm^,Pl>{a}))dJ(a). 



, CeAmb 

No we deduce from Lemma [5.41 that, for all a — (w,j,t) e 1 d x I x K such that 
t < 0, the event {tv(o) C Amb„('P U {a})} is Tt~ measurable. On the other hand, 
for all (x, t) £ 7r(a), the random variable g x j is measurable with respect to o~(V<:t), 
where V<% := V (1 (Z d x Xx] — oo, t[). As a consequence, l(7r(a) C Amb n (V U {a}) 
and t ) £7r ( Q ) 9x,t(V U {«}) are independent. Moreover, one has that 

\(x,t)€7r(a) / (:c,t)G7r(a) 

with 

¥>(*,*) ~ / /(C)d(*(x,t)*M)(C). 

Thus 

(20) E(*(a,?U{«})) =E ( l(7r(a) C Amb„(P U {a})) x ^ ¥>(M) 

\ (x,i)G7r(a) 

Applying again Campbell's theorem, we deduce from ([20| that 



f E(*(a,?U{a}))rfJ(a) = E [ ^] l(7r(a) C Amb„) ^ ^(ar,t) 

\oeP (x,t)eir(a) j 

\CSAmb„ . 

Assuming the conclusion of the proposition to be true for n, we deduce that 
^(a,VU{a}))dJ(a)= f ^(C)rfM*"(C) = / /(C)dAf*(" +1 >(C). 



In view of (fP9l) . this establishes the conclusion of the proposition for n + 1. □ 



Proof of Corollary ] 5. 3[ For n = the result is immediate. For n > 1, 

E /(o^°T C = E*'(« J n 

CSAmb„ o£P 

with 

*'(a,w) := 1(tt(o) C Amb„(w)) ^ fl£ |t (w), 

(x,t)ETr(a) 

Then Campbell's theorem shows that 



E f E fiO-DorA = f E(tt'(a, 



PU{a}))dJ(a). 
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As in the proof of Proposition [521 given a — (w,j, t) eZ^xlxl such that t < 0, 
the event {Tr(a) C Amb n (V U {a})} is J 7 *— measurable while, for all (x,t) S 7r(a), 
the random variable g' x t is measurable with respect to a(P < t). Moreover, 

E{g' xA (PU{a})) = f(x,t)-E(D), 

so that 

E(y'(a,VU{a})) = E I l(?r(a) C Amb„(a;)) ^ f{x,t)-E(D) 
Another application of Campbell's theorem yields that 



/ 



E(V'(a,VU{a}))dJ(a) = E(L») • E ( ^ l(vr(a) C Amb„) f( x ^) 

= e(d)-e\ /(C) 

yceAmb„ 

Proposition 15.21 then yields the conclusion. □ 

Proof of Theorem\^\ Assume that g < 1, and note that, by definition, one has 
J dM(Q — g. We now use Proposition 15.21 with / = 1, and obtain that, for all 
n > 0, E(|Amb n |) < g". As a consequence, 

EflAmbool) <E ^|Amb n | = ^E(|Amb„|) < ^ g" < +oo. 

\n>0 y n>0 n>0 

It is now clear that P(|Amboo| < +oo) = 1. Similarly, applying Corollary 15.31 with 
/ = 1 and D := 1(T = -co) yields that, for all n > 0, 

e ( i^ ^,* = -°°)) =°- 

\(l,t)6Amb„ J 

As a consequence, with probability one, T o T x t > -co for all (x,i) £ Amboo. We 
have thus proved that P(T* > — oo) = 1. The conclusion of the theorem is now a 
consequence of Proposition 15.11 □ 



Proof of Theorem fCl We apply Corollary 15.31 with f(x,t) = exp(Ai) and D := 
exp(AT). As a result, for all n > 0, 



\(i,t)eAmb„ / ^ 

One has 



e xt dM* n (x,t) = (J e xt dM(x,t)j = (A Httime (X)) n . 
Summing over n > 0, we obtain that 

E (e XT ) 



E 



e^°^ X (e AT ) (A ff , time (A))« = . 
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Now by definition of T*, using the fact that A < 0, 



E(exp(AT*)) < E ^ e Kt+To Ttc , t ) 




□ 



Proof of TheoremW\ (sketch). First note that, when g < 1, L* + and L*_ are a.s. 
finite, using an argument similar to the one establishing that T* is a.s. finite in 
the proof of Theorem [A] Moreover, the proof of the estimates on E(exp(AL!j_)) 
and E(exp(— XL*_)) is completely similar to the proof of Theorem [Cl It remains 
to prove that L* indeed bounds the width of T*. This is done by adapting the 
proof of Proposition 15.11 as follows. We work on the event T* > — oo. First note 
that, thanks to the fact that L defines a stopping box, L is T~ L ,L —measurable 
for any random variable such that L < L' . Then observe that Amboo is measurable 
with respect to F~ L ' L , since, for each n > and (x,t) G Amb„, one has that 
L*_ < x — L o r x ,t < x + L o T x _ t < L* + . We now start the induction with the 
case (Amboo | = 1< Then the values of the e(a, T, £) for a G H are completely 
determined by H itself, and H is measurable with respect to F~ L,L , so we are 
done. If |Amboo| = n + 1, we apply the induction hypothesis to every |Amboo °T x ,t\ 
such that x — z + y for some a = (z, i,t) G Ambi and y G A,-, then deduce that 
(0](0) has the required measurability properties. □ 



We start with a proposition showing that (T, H) is a CFTP time with ambiguities 
for the perturbed dynamics if Xoo is a finite set P— a.s. 

Proposition 6.1. IfF^X^] < +oo) = 1, then (T, H) is a CFTP time with ambi- 



Proof. First note that the stopping property of I m is a direct consequence of the 
way the process is constructed. We now work on the event that \3Zoo\ is finite. 
Consider £ G S z , and let Q be the element of fl u obtained from V by replacing 
any a — (x, i,s) G V such that s >T and i G T 9 by (x, t e ( a .T,£) , s), and suppressing 
any a — (x,i,s) G V such that s < T and i G X p . Let denote the random 

variable defined on VL U by 



where denotes the stochastic flow defined by V w on 17". From the definition of 
the dynamics, we see that 



Now define a path cq, . . . , c m in T as follows. Start with cq := r. Then assume that 
Co, . . . , Cfe have been defined. If Ck has no son in T, the path ends at Ck, so that 
m := k. If Cfe has a single son in T, then Ck+i is defined to be this single son. Finally, 
if Ck has several sons in T, Ck+i is defined to be the son associated with the value 
v := e((x Ck , i Ck , t Ck ), T, £). By definition of the exploration processes X" and X, one 



has that X^(Q) = ± Cm . Then by definition T U (Q) = mf{t; (x,i,t) G X^(Q)}, so 
that the identity X^(S) = X Cm implies that T U (Q) > T. From the fact that T u 



6. Proof of Theorem IbI 



guities. 



h t (0 := [($") t °-(0](0), 



(21) 



[*t (01(0) = [MOKS). 
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is associated with the exploration process specified by 6, there exists a measurable 
map G defined on the set such that, on the event {T u > —00} = {|X^| < +00}, 

one has, for every \ G S zd , hr^(x) = G(X!^). As a consequence, 

(22) [M£)KQ) = G(£cJ. 

It is now immediate from (f^Tj) and ((^)) that if e(a,£i,T) = e(a,^2,r) for every 
a £ H, then = [^t~(£2)](0), since both values of X Cm obtained starting 

from £1 or £2 are identical. □ 

The next two propositions are the key first-moment estimates needed to control 

Define a kernel K on flj as follows. If 6(X) = 0, then K(X, •) = 5 X (-)- If 0(X) ^ 0, 
let s :— inf{i; (x, i, t) G X}, and let 

dAT(X,XU{(x,i,i)}) =nexp(|0(X)|r„(t-s))l(t < s)l(x G 9(X))dJ u (x,i,t), 

with r u := X^ez u r i- From the definition, one has the following. 

Proposition 6.2. The sequence (Xg)e>o is a Markov chain on with initial state 
and transition kernel K . 

Now define a kernel K on the set of elements lellj such that 9(X) ^ as follows. 
Let s := inf{t; (x,i,t) G A}, and let 

dK(X, X U {(a;, i, i)}) = f; exp(|0(X)|t(f - s))l(i < s)l(x G 6(X))dJ u {x,i,t), 

with, for i G X u , fj := rj + ^2 jeXP X^g/^a,) r jH L v = i), and r := J2 3 ei r r Denne 
also the kernel L by 

dL(X, {x,i,t)) = nex.p(\0(X)\t(s - t))l(t < s)l(x G 6{X))dJ{x,i,t). 

For £ > 0, let IV (resp. T^) denote the set of vertices at distance I from the root in 
T (resp. T'); let also 

T t := {(X , . . . , Xg) G (SIJ)' +1 ; 9(X ) + 0, . . . , ^ 0} , 

A, := {(X , . . . , Xg) G (fi]f) <+1 ; 0pf o ) ^ 0, . . . , ^ 0} . 

Proposition 6.3. For every i > 0, and any measurable map F : — > K +; one 

/las </ie identity 

Aor every £ > 0, ana 1 any measurable map / : Z d x I x K — > R + , one /ias i/ie 
ideniiiy 



E ]T /(x a ,z a ,i )] = /. 



f(a)d5® (X )dK(X , Xi) ■ ■ ■ dK{X t -i,Xi)dL{X h a), 
where the integral is over (Xq, . . . , Xi, a) G Ag x (Z rf xIxR). 
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Proof. The proof is similar to that of Proposition 15.21 We start with the first 
identity, whose proof is by induction. For I = 0, 1, the result is a direct consequence 
of the definition. Now for £ > 1, write 

(23) n*a) = ^i(a,n 

where 

Z x (ev):= l(a = (z c ,*c,*c)) 51 •), 

with b <— c meaning that 6 is a son of c in T, and with 
Z 2 (X r ):=l{jGl u )F(XU{(yJ,s)}) + l(jGl p ) £ F(X U {(„, i v , a)}), 

with (y,j,s) := n(9(X),t), and i := inf{^; {w,k,v) e X}. (We write Zi(a, •) and 
•) to make the dependence on uj explicit.) 

By Campbell's theorem (Theorem |F|) . one has that 

E = f^{Zi(a,PU{a}))dJ(a). 

Now given a = (x,i,t), define the random finite counting measur 

M:= l(« = (^,ie,tc)) E lW*6)^0)«5 i6 , 

so that 

Zi(a,-) - / Z 2 (X, -)dM(X, •). 



Note that .M is J^— measurable, while, for any X G 51^ such that #(X) 7^ and 
t := infjV; (w,k,v) G X}, Z 2 {X, •) is measurable with respect to o~(V<:t) and 
satisfies E(Z 2 (X, •)) = KF(X). We deduce that 



E(Zi(a,^U{a})) =E (/ #F(X)dM(X, V U {a}) J = E(Z 3 (a, P U {a})), 
where 

Z 3 (a,-):= l(a=(xc,ie,t c )) J2 ^F{± h ). 

c€T e _ 1 beTi, b-^c 

Using Campbell's theorem again shows that 

f E(Z 3 (a, V U {a}))dJ(a) = E ( £ KF(X b ) j . 

This computation allows induction over £ to be used to prove the desired identity 
for all £ > 1. The second identity of the Proposition can be deduced from the first 
one, using an argument similar to the derivation of Corollary 15.31 from Proposition 

EI □ 



^We equip the set of finite counting measures on with the a— algebra generated by all the 
maps of the form Y h> Y(B), where B belongs to the cr— algebra defined on O^. 



22 



JEAN BERARD, DIDIER PIAU 



Combined with Propositions 16.21 and 16. 3[ the following remark is the key to ob- 
taining estimates on X from the properties of X u . Consider 1 6 Uj such that 
6(X) ^ 0. From the definition of e given in ([§]) and the fact that x u < r, one has 
that 

(24) dK{X, X U {{x, i, t)}) < (1 + e)dK(X, X U {(x, i, t)}). 

We can now prove the various estimates that are needed in the proof of Theorem 

m 

Lemma 6.4. For every £>1, 

E(\T e \)<(l + e) e P u (\X^\>£). 

Proof. By Proposition 16. 31 one has that 

E(|T,|) = / d6f l (X )dK(Xo,X 1 )---dK(X^ 1 ,X e ) 1 
J(x ,...,x e )evt 

From (|24| . one deduces that 



H\T e \) < {1 + eY / dS ti {X )dK{X ,X 1 )---dK(X i _ 1 ,X i ). 
J(x ,...,x t )er e 



But by Proposition 16. 2\ 

d6 $ (X )dK{X ,X 1 )---dK(X e _ 1 ,X i )=F u (\Xl > \>£), 



i(x ,...,x e )er e 
so that 

E(|T,|)<(l + e) f P"(|X^|>£). 

□ 

Lemma 6.5. For every £ > 0, 

El ^|AJl(i Q 6^)J < K (l + eYF u (\X^\>£). 

\a€T e J 

Proof. Let f(x, i, t) := \Ai\l(i £ IP). Note that, given X G such that B{X) ^ 0, 
one has J^ x . t ^^ dKlKK f{x,i,t)dL(X, (x,i,t)) — k. Using Proposition ^. 31 and (l24l) . 
one obtains that 

El \A la \l(iatl p )) <k(1 + £) £ / 2 , 

where 

h= f d6 $ (X )dK(X ,X 1 )---dK(X i _ 1 ,Xe)=F u (\X^ o \>£). 
J(x ,-,x e )eA e 

□ 

Lemma 6.6. For every £ > 0, and A < t, 

E I J2 ex P (Ai a ) ) < ^(1 + e) £ E"(exp(A^_ 1 )l(|X^| > £)), 
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with the convention t-% := (remember that (x n ,i n ,t n ) n >o corresponds to the 
exploration process X u ). 

Lemma 6.7. For every £ > 0, and A < t, 



E \Ai\exp(Xt a )l(i a eTP) < ^-^(l + e)*E u (exp(At^i)l(|X£,| >£)), 

with the convention f_i := 0. 

Proof. We prove Lemma 16.71 the proof of Lemma 16.61 being quite similar. Let 
f(x,i,t) := L4i|exp(A«)l(i e l p ). Consider X € flj such that 6{X) ^ 0, and 
let s ;= mi{t; (x,i,t) 6 X}. One has j {x i t)eJjdxIyM f(x,i,t)dL(X,(x,i,t)) = 

Kx\9(X)\(t\6(X)\ - A)" 1 exp(As) < Kt(t — X) 1 exp(As), since \9(X)\ > 1. Using 
Proposition 16. 31 and (f24|). we deduce that 



E ]T exp(At a )l(z a £P) < kx(x - A) _1 (l + e)%, 

with 

h = I ex V (Xs^ 1 )dS (X o )dK(X o ,X 1 )---dK(X i . 1 ,X i ) 
J{Xo,- ,x e )eA e 

and S£-i :— inf{£; {x,i,t) G Xi}, with the convention s^_i := 0. Then note that 
7 3 = E u (cxp(A* < _i)l(|^| >£)). □ 

Lemma 6.8. For every £ > 1, A E K and 1 < q < d, 



e (E E <o(i+ e )^«( 1 ^ r £ 

\iei p zeAi ) 



e Ax «l(|*« \>t) I , 



where 
(25) 



We are now ready to prove Theorem [Bj From now on, we assume that there exists 
fj, > such that E u (e^' 3 ^ |) < +oo. As a consequence, there exists a finite constant 
C such that, for all t > I, 

(26) P"(|X^| >f) <Ccxp(-^). 

We first prove that, as soon as e < fi, E(|j£oo|) < +oo. Indeed, one has that 
|Xoo I < |T|, so that 

(27) E(|i 00 |)<E(|T|)=El^|T,|] = $>(|T,|). 

By LemmaEH for all I > 1, one has E(|T £ |) < (1 + e) f F"(|X^| > I). Combining 
(|26|) and {27}, we see that E(|Xoo|) < +oo when e < //. By Proposition 16.11 this 
proves that the pair (T, H) defined by (TTT|) and (TT2|) is indeed a CFTP time with 
ambiguities. 
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We now prove that, for small enough e and K, the pair (T, H) satisfies < 1. Using 
the definition, then Lemma 16.51 , we have that 



= E E £i^.i 1 & e2P ) <Kj2(i+zfv u (\x u 00 \>t)- 

e>o \a£T e J e>o 

From (|26l) . we see that g < 1 for all e < ijl and small enough k. 

We now prove that E(exp(AT)) < +oo for all small enough e and A. We start with 
the observation that 

E(exp(AT)) < E ( ^ exp(At a ) J = E E E e M^ a ) 



From Lemma Win we deduce that, for all A < r, 

(28) E(exp(AT)) < ^(1 + e)^E u (exp(At,_ 1 )l(|X^| > t)). 

By Schwarz's inequality, 

E"(exp(At,_ 1 )l(|^ | > £)) < (E it (ex P (2A^_ 1 )) 1/2 P"(|X^| > i) 1 ' 2 . 

Bounding above by the sum of I independent exponential random variables 
with parameter r u on one hand, and using (|26l) on the other hand, one obtains that, 
when A < r u /2, 

(29) E"(exp(A^)l(|X^ > ^)) < C 1 / 2 exp(-^ + l)/2). 

Combining (|2"51) and (|29l) . we have that E(exp(AT)) < +oo for all small enough e 
and A. 

We now prove that A# time (A) < 1 for all small enough e, k, A. Using Lemma 16.71 
we obtain that 

(30) A ff , time (A) < V(l + e/E"(cxp(A^)l(|X l 4| > I)). 

Using again (|29|) . one concludes that A#,ti„ le (A) < 1 for all small enough e, k, A. 

Now let R denote the depth of T, and define L :— (3 U (R) (remember that (3 is 
defined in ©). By definition of the exploration process with locking of ambigu- 
ities, one checks that L defines a stopping box and that H is measurable with 
respect to F~ L > L . Now, <(2T]) and (J22|) show that (0) satisfies the required 

measurability properties. 

Using the obvious inequality R < |T|, Lemma [6.41 shows that A^(A) is finite for 
small enough e, A. Finally, Lemma |6 . 81 shows that A# )Space (A, q) < 1 for all q, when 
e, k, A are small enough. 
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